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To clarify the effect of deep rooting on grain yield in rice ( Otyza sativa L.) in an irrigated paddy field with or 
without fertilizer, we used the shallow- rooting IR64 and the deep-rooting Drol-NIL (a near-isogenic line 
homozygous for the Kinandang Patong allele of DEEPER ROOTING 1 (DROl) in the IR64 genetic 
background). Although total root length was similar in both lines, more roots were distributed within the 
lower soil layer of the paddy field in Drol-NIL than in IR64, irrespective of fertilizer treatment. At maturity, 
Drol-NIL showed approximately 10% higher grain yield than IR64, irrespective of fertilizer treatment. 
Higher grain yield of Drol-NIL was mainly due to the increased 1000-kernel weight and increased 
percentage of ripened grains, which resulted in a higher harvest index. After heading, the uptake of nitrogen 
from soil and leaf nitrogen concentration were higher in Drol-NIL than in IR64. At the mid-grain-filling 
stage, Drol-NIL maintained higher cytokinin fluxes from roots to shoots than IR64. These results suggest 
that deep rooting by DROl enhances nitrogen uptake and cytokinin fluxes at late stages, resulting in better 
grain filling in Drol-NIL in a paddy field in this study. 



Plant roots are important organs for the uptake of water and nutrients from the surrounding soil. The 
physiological and morphological traits of roots affect shoot growth. In rice (Oryza sativa L.), a well- 
developed root system increases biomass and yield in different treatments and cultivars in paddy fields 1 ". 
Higher uptake of nitrogen and water during the grain filling stage allows the plant to attain higher dry matter 
production by maintaining a higher photosynthetic rate in paddy fields 1 ' 5-7 " 9 . High flux of cytokinins, which are 
synthesized in the roots and transported through the xylem to the shoot 1014 , contributes to the maintenance of 
high photosynthetic rate owing to increasing leaf nitrogen and ribulose- 1,5-bisphosphate carboxylase/oxygenase 
(Rubisco) which is a key enzyme in photosynthesis in rice 6,9-15,16 . High-yielding rice plants often have a larger 
volume of downward-growing nodal roots, which increases the proportion of roots in the lower soil layer in paddy 
fields 17-20 . However, the importance of these root ideo types has been suggested on the basis of comparison among 
cultivars with different root and shoot morphologies, and it is not clear whether the deep root system contributes 
to high yield and dry matter production in paddy fields. To accurately evaluate the effect of deep rooting on rice 
grain yield, homogeneous genetic materials such as near-isogenic lines (NILs) with differences only in the root 
growth angle are needed. 

Recently, Uga et al. 21 cloned and characterized a quantitative trait locus (QTL) for root growth angle in rice, 
DEEPER ROOTING 1 (DROl). DROl functions downstream of the auxin signaling pathway and controls the 
gravitropic curvature in rice roots. The authors also developed Drol-NIL, carrying a functional allele of DROl 
derived from the deep-rooting cultivar Kinandang Patong (a traditional tropical japonica upland cultivar of 
Philippine origin) in the genetic background of the shallow- rooting parent variety IR64 (a modern lowland indica 
cultivar, widely grown in South and Southeast Asia), which has a non-functional allele of DROl. Under upland 
conditions, Drol-NIL showed deeper rooting than IR64, although the two lines did not show any marked 
differences in other root and shoot traits 21,22 . Under upland conditions with drought stress, Drol-NIL maintained 
better grain filling and thus showed higher yield than IR64 because the plants avoided drought stress by collecting 
water from deep soil. This clearly demonstrated that deep rooting is a beneficial trait for stable rice production 
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Table 1 | Mean total root length and root length density in different 
soil layers in IR64 and Drol -NIL lines as affected by line and fert- 
ilizer treatment 



Line (L) 

IR64 
Drol-NIL 



ANOVA 

L 

F 

LxF 



Root length density (cm cm 3 

Fertilizer Total root length 

treatment (F) (cm core-') 0-5° 5-10° 10-20° 



No treatment 
Treatment 



3199 
3268 
3160 
3307 

ns 
ns 
ns 



20.3 
16.1 
18.0 
18.3 



ns 
ns 



6.5 
10.0 
7.8 



ns 
ns 



2.9 
3.6 
3.2 
3.3 



ns 
ns 



"Distance from the soil surface (cm). 

Root length density was sampled near hills between rows at 20 days after heading. Only >0. 1-mm- 
thick roots were measured. 

"Significant at P < 0.01; ns, not significant at P< 0.05. P values are based on ANOVA. 



under drought conditions. Thus, we assumed that Drol-NIL might 
be a useful material to clarify the effect of deep roots in paddy fields. 

Lodging also affects grain yield and dry matter production 23 . 
Although previous studies examined the relationship between deep 
rooting and lodging resistance, they compared different cultivars or 
the effects of root pruning 24 27 . Such comparisons cannot exclude the 
effects of genotype differences or root injury, respectively, on the 
observed effects. 

The objective of this study was to examine the effect of deep root- 
ing on grain yield in paddy fields by using the deep-rooting Drol- 
NIL and the shallow-rooting IR64. We evaluated not only yield and 
its components, but also physiological traits such as nitrogen uptake 
from roots and cytokinin fluxes in xylem exudates and leaves. We 
also investigated the effect of deep-rooting Drol-NIL on lodging 
resistance. 

Results 

Root length density. There was no significant difference in total root 
length between IR64 and Drol-NIL, irrespective of fertilizer 
treatment (Table 1). In both lines, root length density decreased 
with increasing depth. Root length density differed significantly 
between the lines when measured in specific soil layers, although 
there was no interaction between the effects of the line and 
fertilizer. Root length density of Drol-NIL was 0.79 times that of 
IR64 in the upper layer, 1.54 times in the middle layer, and 1.24 times 
in the lower layer, irrespective of fertilizer treatment. 

Dry matter production, yield and yield components. The dry weight 
of the above-ground parts was similar in both lines at heading, but was 
higher in Drol-NIL at maturity, irrespective of fertilizer treatment 
(Table 2). Grain yield and HI were also higher in Drol-NIL than in 
IR64, irrespective of fertilizer treatment. The number of panicles per 
m 2 , the number of spikelets per panicle and the number of total 
spikelets per m 2 were similarly increased by fertilizer treatment in 
both lines (Table 3). In contrast, the 1000-kernel weight and 
percentage of filled spikelets were higher in Drol-NIL than in IR64, 
irrespective of fertilizer treatment. 

Nitrogen concentration in leaves and nitrogen content of plants. 

At heading, nitrogen concentration in flag leaves and total nitrogen 
content of plants were similar in the two lines (Tables 4 and 5). On 
the other hand, at 20 days after heading, nitrogen concentration in 
flag leaves was higher in Drol-NIL than in IR64, irrespective of 
fertilizer treatment (Table 4). At maturity, total nitrogen content 
of plants was also higher in Drol-NIL than in IR64, irrespective of 
fertilizer treatment (Table 5). As a result, nitrogen accumulated 



Table 2 | Mean dry matter accumulation, grain yield, and harvest 
index in IR64 and Drol NIL lines as affected by line and fertilizer 
treatment 



Line (L) 


Fertilizer 
treatment (F) 


DWH° 
(t ha" 1 ) 


DWM° 
(tha- 1 ) 


Yield b 
(t ha" 1 ) 




IR64 




10.92 


15.47 


6.15 


0.43 


Drol-NIL 




1 1.24 


16.27 


6.80 


0.45 




No treatment 


8.86 


12.21 


4.71 


0.43 




Treatment 


13.30 


19.52 


8.24 


0.45 


ANOVA 












L 




ns 


* 


** 


** 


F 




** 


** 


** 


ns 


LxF 




ns 


ns 


ns 


ns 



°DWH, dry weight at heading; DWM, dry weight at maturity. 
b Yield is given for unhulled rice with water content of 1 5%. 

c Harvest index (HI) was calculated by dividing the dry weight of panicle by the dry weight of the 
above<3round parts. 

♦Significant at P < 0.05; ** significant at P < 0.0 1 ; ns, not significant at P < 0.05. Pvalues are 
based on ANOVA. 



during grain filling was higher in Drol-NIL than in IR64 with or 
without fertilizer treatment, although nitrogen concentration in 
leaves and nitrogen content of plants were increased by fertilizer 
treatment in both lines (Tables 4 and 5). 

Xylem exudation rate, cytokinin flux in xylem exudates and 
cytokinin concentration in flag leaves. At 14 days after heading, 
the rate of xylem exudates and the flux of cytokinins (determined as 
frarcs-zeatin, ris-zeatin, frarcs-zeatin riboside, N 6 (A 2 -isopentenyl) 
adenine, N 6 (A 2 -isopentenyl)adenosine, dihydrozeatin, and dihydro- 
zeatin riboside combined) in xylem exudates, as well as concentra- 
tions of these cytokinins in flag leaves, were higher in Drol-NIL than 
in IR64, irrespective of fertilizer treatment, although both parameters 
were increased by fertilizer treatment in both lines (Table 6). 

Lodging resistance. At 14 days after heading, pushing resistance was 
higher in Drol-NIL than in IR64 (Table 7). The difference between 
the two lines was no longer significant at 35 days after heading. 

Discussion 

The roots of Drol-NIL were distributed in the lower soil layer than 
those of IR64 in a paddy field (Table 1 , this study) , similar to previous 
observations in upland fields 21 . Recently, we have demonstrated that 
Drol-NIL produces more filled grains than IR64 under drought 
stress conditions, but not under non-drought conditions 21 . These 
results were obtained in upland fields. Thus, we wondered whether 
yield performance differs between shallow- and deep-rooting vari- 
eties in paddy fields. 

In a paddy field, grain yield of the Drol-NIL was higher than that 
of IR64, irrespective of fertilizer treatment (Table 2). Under drought 
conditions, deeper rooting caused by DROl led to a lower percentage 
of seed sterility than in shallow-rooted IR64 and thus resulted in 
higher yield, because plants avoided drought stress by increasing 
water capture from deep soil 21 . In a paddy field, higher grain yield 
of the Drol-NIL was mainly due to the higher 1000-kernel weight 
and percentage of ripened grains, resulting in better grain filling and 
higher HI (Tables 2 and 3). This suggests that the Drol-NIL may 
have greater ability to supply carbohydrates to the panicle after 
heading. 

Plant carbohydrates are produced by photosynthesis. Nitrogen is 
essential for the production of photosynthetic enzymes, and there is a 
positive relationship between leaf nitrogen content and the rate of 
photosynthesis 6,9,28 . Recently, Takai et al. 29 provided genetic evidence 
suggesting that photosynthetic rate is positively related to leaf nitro- 
gen content through the cloning and characterization of GPS 
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Table 3 Mean yield components in IR64 and Drol -NIL lines as affected by line and fertilizer treatment 






Panicles 


Spikelets per panicle 


Total spikelets 


1 000-kernel weight" 


Filled spikelets b 


Line (L) 


Fertilizer treatment (F) 


(no. m" 2 ) 


(no.) 


(X10 3 m" 2 ) 


(g) 


[%) 


IR64 




280 


106 


29.9 


23.1 


0.68 


Drol-NIL 




282 


104 


29.5 


23.6 


0.73 




No treatment 


261 


91 


23.6 


23.4 


0.69 




Treatment 


301 


119 


35.7 


23.3 


0.72 


ANOVA 














L 




ns 


ns 


ns 


** 


** 


F 




** 


** 


** 


ns 


ns 


LxF 




ns 


ns 


ns 


ns 


ns 


-1,000-kernel we 


ight is given for water content of 1 5%. 












Percentage of fil 


ed spikelets refers to brown rice with a c 


rain thickness of > 


.6 mm, 








"Significant at P< 0.01; ns, not significant at P< 0.05 


P values are based 


on ANOVA. 









(GREEN FOR PHOTOSYNTHESIS), a rice QTL that regulates car- 
boxylation efficiency and thus controls the photosynthetic rate. In 
this study, nitrogen concentration in flag leaves (Table 4) and the 
other leaves (data not shown) after heading were higher in Dro 1 -NIL 
than in IR64, irrespective of fertilizer treatment, which may conduce 
better grain filling in Drol-NIL. We hypothesize that increase in 
nitrogen uptake by DROl enhances photosynthetic rate, resulting 
in higher grain filling ability. However, we did not investigate the 
photosynthesis of leaves in both lines of this study. Therefore, the 
relationship between photosynthesis and grain filling remains to be 
investigated. 

Nitrogen is supplied to rice from both fertilizer and mineralized 
soil nitrogen 30 . Because nitrogen in the upper soil layer is often insuf- 
ficient at the late growth stage, efficient uptake of nitrogen from the 
lower layer is important for vigorous growth at this stage 31 . A sig- 
nificant difference in plant dry weight and nitrogen uptake between 
IR64 and Drol-NIL at maturity, but not at heading, suggests that 
higher root density of Dro 1 -NIL in the lower soil layer in comparison 
with IR64 was probably advantageous in capturing soil nitrogen 
from these layers during the grain filling period. However, the rate 
of nitrogen uptake from the lower soil layer still needs to be accur- 
ately measured to clarify how nitrogen absorbed from the lower soil 
layer in rice carrying DROl assists increasing grain yield. We 
observed several roots of Drol-NIL penetrated through the plow 
pan layer (data not shown), although we did not quantify roots 
and nitrogen existed below 20 cm in the soil layer. We also have to 
elucidate effect of deep roots below 20 cm of soil layer on nitrogen 
uptake. 

Cytokinins are involved in the control of leaf senescence 32-33 . The 
delay in leaf senescence is considered to be important for increased 



Table 4 | Mean nitrogen concentration in flag leaves in IR64 and 
Drol NIL lines as affected by line and fertilizer treatment 



Line (L) 



Fertilizer 
treatment (F) 



Nitrogen concentration (g g ') 



Headii 



20 DAH° 



IR64 
Drol-NIL 



ANOVA 

L 

F 

LxF 



No treatment 
Treatment 



0.024 
0.024 
0.022 
0.026 



0.017 
0.019 
0.015 
0.020 



ns 



Q DAH, days after heading. 

•Significant at P< 0.05; " significant at P < 0.01; ns, not significant at P< 0.05. P values are 
based on ANOVA. 



dry matter and grain production in paddy rice 1 ' 5 7,9,34 . Higher trans- 
port of cytokinins from roots to shoots during grain filling is also 
associated with the maintenance of higher leaf nitrogen and Rubisco, 
and thus higher photosynthetic rate 6,9,15,16 . In this study, both cyto- 
kinin flux in exudates and cytokinin concentration in flag leaves were 
higher in the Drol-NIL than in IR64, irrespective of fertilizer treat- 
ment (Table 6). These results indicate that DROl plays an important 
role in elevating the cytokinin level. We suppose that this higher 
cytokinin level of Drol-NIL was associated with maintenance of 
higher leaf nitrogen and delaying leaf senescence. How is Drol- 
NIL able to maintain higher cytokinin levels during grain filling? 
The roots distributed in the lower layer were white and fresh in 
comparison with the roots distributed in the upper soil layer, which 
implies a non-reducing environment around the plow pan layer. We 
also observed that roots of Drol-NIL penetrated through the plow 
pan layer were white and fresh (data not shown). Therefore, we 
hypothesize that elevated cytokinin synthesis in Drol-NIL after 
heading was due to its vigorous roots in the lower soil layer and 
the layer beneath the plow pan. The mechanisms of high cytokinin 
flux in rice plants with roots distributed in the lower layer and its 
effects on photosynthesis and yield performance remain to be 
clarified. 

Lodging reduces grain yield and dry matter production in paddy 
rice 35 . In this study, Drol-NIL showed higher pushing resistance as 
an index of lodging resistance than IR64 before the soil hardness 
increased as a result of the drainage of residual water in the paddy 
field, although lodging no longer occurred after heading in either line 
(Table 7). This suggests that roots distributed in the lower layer 
contributed to the increased root lodging resistance in wet soil. 
Because the position of the shoot base of rice plants directly sown 
in the soil is shallower than that of transplanted rice plants, lodging 
easily occurs after heading in direct-sowing cultivation 27 , resulting in 
a reduction in grain yield and dry matter production 23 . The roots 
penetrating into the lower soil layer with a higher bulk soil density 
have been proposed to be an important determinant of lodging 
resistance in rice in direct sowing cultivation 25 27 . We have now 
demonstrated that introduction of the deep rooting trait into shal- 
low-rooting cultivars can form the basis of a breeding strategy for 
lodging resistance. 

Under upland conditions, the maximum root depth of the Drol- 
NIL was more than 40 cm while that of IR64 was around 20 cm 21 . 
The plow pan was 20 cm below the soil surface in our paddy field. 
The root growth of Drol-NIL might be affected by the plow pan but 
some roots penetrated through the plow pan layer. Since we did not 
measure a penetration ability of roots in both lines, whether DROl 
improved the penetration ability of roots is unclear. It is necessary to 
further investigate the effect of DROl on penetration ability to deter- 
mine if roots that penetrate through the plow pan have some benefits 
with respect to nitrogen uptake and cytokinin fluxes. 
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Table 5 Mean nitrogen content of whole plants in IR64and Drol - 
NIL lines as affected by line and fertilizer treatment 

Total nitrogen content (g m~ 2 ) 
Fertilizer Accumulated 

Line (L) treatment (F) Heading Maturity nitrogen" (g m" 2 ) 


IR64 


8.6 


11.5 2.8 


Drol -NIL 


9.1 


12.3 3.2 


No treatment 


6.0 


8.1 2.1 


Treatment 


1 1.7 


15.7 3.9 


ANOVA 






L 


ns 


** * 


F 


** 


## ## 


LxF 


ns 


ns ns 


"Accumulated nitrogen was calcL 
maturity. 

•Significant at P < 0.05; ** sign 
based on ANOVA. 


lated as the difference between the total content at heading and 
ficant at P < 0.01 ; ns, not significant at P< 0.05. P values are 



Our results suggest that deep rooting caused by the functional 
allele of DROl improves grain yield in a paddy field by enhancing 
the higher 1000-kernel weight and percentage of ripened grains. 
DROl predominantly influences root growth angle rather than other 
root morphological traits such as root length and number 21 . 
Physiological and yield traits improved in the Drol -NIL should be 
attributable to change of vertical root distribution by DROl , although 
we have to elucidate relationships between higher nitrogen uptake 
and cytokinin fluxes by DROl and grain yield in the Drol -NIL. 
Previous studies have reported QTLs associated with photosynthesis 
and grain yield near the DROl region on chromosome 9 3637 . Multiple 
QTLs for net photosynthesis rate were found around SSR marker 
RM410 at flowering and grain filling under both drought and well- 
watered field conditions 36 . Dixit et al. 37 detected promising QTL, 
qDTY 91B , for grain yield under paddy fields with drought stress in 
the interval between SSR markers RM24350 (15.4 Mb) to RM24390 
(15.9 Mb), whereas this QTL did not show positive effect on grain 
yield under non-drought conditions. They supposed that qDTY 9 !B 
was due to pleiotropic effect of DROl. The IR64 allele of DROl was 
identified in several IR64 progenies showing shallow rooting but not 
in the ancestors of IR64 21 . Introducing the Kinandang Patong allele 
of DROl into these modern varieties might be improved their grain 
yield in the paddy fields, although potential effects of DROl on grain 
yield in other genetic backgrounds remain to be clarified. Thus, 



Table 7 | Pushing resistance in IR64 and Drol NIL lines 

Pushing resistance (kg) 



Line(L) 1 4 DAH" 35 DAH° 



IR64 


1.32 


1.30 


Drol -NIL 


1.45 


1.35 


ANOVA 






L 


* 


ns 



a DAH, days after heading. Ponding water was drained starting from 25 days after heading. 
*SignificantatP< 0.05; ns, not significant at P< 0.05. Pvalues are based on ANOVA. 



further analysis using the other genetic backgrounds will be needed 
to understand relationship between deep rooting by DROl and its 
pleiotropic effect. 

In the developing countries lacking agricultural infrastructure, 
farmers may not be able to use sufficient amounts of nitrogen ferti- 
lizers in paddy fields. DROl may be useful in these areas because 
nitrogen uptake from the lower soil layer due to deep rooting would 
contribute to increasing grain yield. However, these effects by DROl 
would depend on the soil fertility and the method of fertilizer 
application in each paddy field. Multi-location trials for the Drol- 
NIL will be needed to conclude whether deep rooting by DROl shows 
universally effective for grain yield in the paddy fields because our 
results were obtained from the one location with relatively small plot 
size. 

Methods 

Materials and plant cultivation. IR64 and Drol -NIL were cultivated in a paddy field 
at the NICS in Tsukubamirai (lat 36°02'N, long 140°04'E), Ibaraki, Japan, in 2010 
and 2011. Drol-NIL was developed from a cross between IR64 and Kinandang 
Patong (as a donor of the functional allele of DROl) by repeated backcrossing with 
IR64 and marker-assisted selection to eliminate non-target DROl regions 21 . The 
Drol-NIL had only homozygous for the Kinandang Patong allele between 15.9 Mb 
(SSR marker RM24386) to 18.8 Mb (SSR marker RM242) on chromosome 9, 
indicating that more than 99% of genome in the Drol-NIL was derived from IR64. 
Plants were grown at a density of 22.2 hills (15 cm X 30 cm) per m 2 , with one plant 
per hill. Germinated seeds were sown in nursery boxes in April and grown in a vinyl 
shelter. In May, when the fourth leaves were expanding, seedlings were transplanted 
into an irrigated paddy field after puddling. The experimental plants were laid out in a 
split-plot arrangement with or without chemical fertilizer, with three replications 
(9 m 2 in 2010 and 28 m 2 in 2011 for each replicate, with fertilizer treatment as the 
main plot and line as the sub-plot). For the fertilizer plot, chemical fertilizers were 
applied as a basal dressing at 9 g m" 2 each of P 2 O s and K 2 0, and 12 g m~ 2 N (6 g 
LP40 and 6 g LP100); LP40 and LP100 are controlled release fertilizers releasing 80% 
of the total nitrogen content at a uniform rate up to 40 and 100 days after application, 
respectively, at 20-30°C. The soil type was the Gleyic Fluvisols, and the soil contained 
5.0 mg 100 g~' available nitrogen in the upper layer (0-10 cm) and 4.4 mg 100 g -1 in 
the lower layer (10-20 cm), as determined by hot water extraction. The plow layer 
was approximately 20 cm from the soil surface, and the lower layer formed a plow 
pan, which was approximately 10 cm thick. Daily mean, maximum and minimum air 
temperatures during growth period were 23.4 'C, 30.4°C and 13.3 'C in 2010, and 
23.0°C, 29.7°C and 13.6°C in 2011, respectively. The heading dates (defined as the 
dates when 50% of plants were at heading) of both lines were August 16, 2010 and 
August 15, 2011 with fertilization, and August 14, 2010 and August 13, 2011 without 
fertilization. Plants were harvested on September 26, 2010 and on September 25, 
2011. 

Measurements of root length density. To evaluate root length density of rice plants 
in paddy field, we used core sampler method, because the core sampler can minimize 
destruction of the soil structure of paddy field compared to other methods such as 
monolith sampler. Four hills with the normal number of panicles were collected from 
each replicate at 20 days after heading in 2010 and 201 1. A soil core (50 mm in 
diameter and approximately 200 mm in depth from the soil surface to the plow pan) 
was collected from position between the rows by using a core sampler (DIK-110C, 
Daiki Rika Kogyo Co., Saitama, Japan). The soil cores were divided into three 
segments corresponding to the upper (0-5 cm), middle (5-10 cm) and lower (10- 
20 cm) layers. The segments were washed carefully to separate roots from soil. The 
root length was determined by using a root scanner (WinRHIZO, Regent 
Instruments, Inc., Quebec, Canada). The root length density in each segment was 
calculated as root length divided by volume of each segment. 

Measurements of yield, yield components and plant dry weight. Plants were 
harvested from an area of 1.1 m 2 in 2010 and 2.6 m 2 in 201 1 for each replicate. Fully 



Table 6 | Exudation rate and mean cytokinin flux in exudates, cyto- 
kinin concentration in flag leaves at 1 4 days after heading in IR64 
and Drol -NIL lines as affected by line and fertilizer treatment 

Cytokinin 

Exudation Cytokinin flux in concentration in 
rate exudates 0 flag leaves 

Fertilizer 

Line(L) treatment (F) (ghill^h" 1 ) (pmol hill" 1 h" 1 ) (pmol g" 1 f.w.) 



IR64 




3.70 


13.81 


2.89 


Drol-NIL 




4.00 


16.27 


3.48 




No treatment 


3.47 


11.90 


3.05 




Treatment 


4.17 


17.86 


3.32 


ANOVA 










L 




* 


* 


* 


F 




* 




* 


LXF 




ns 


ns 


ns 



"The flux of cytokinins (frans-zeatin, c/s-zeatin, rrans-zeatin riboside, N 6 - (A 2 -isopentenyl)adenine, 
N 6 (A 2 -isopentenyl) adenosine, dihydrozeatin, and dihydrozeatin riboside) was calculated from 
the exudation rate and cytokinin concentration. 

'Significant at P< 0.05; ** significant atP< 0.01; ns, not significant at P< 0.05. Pvalues are 
based on ANOVA. 
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ripened grains were selected by sieving through a 1.6-mm mesh, and the moisture 
content was adjusted to 15%. The harvest index (HI) was calculated as the panicle dry 
weight divided by the total plant dry weight. To determine plant dry weight, three (in 
2010) or eight (in 201 1) hills were selected for each replicate. Each plant was separated 
into leaves, culm and leaf sheaths, and panicles. Samples were dried in an oven at 80 L 'C 
for 5 days and weighed. 

Measurement of nitrogen in flag leaves and whole plants. Measurements were 
performed in 201 1. Flag leaves were collected from 8 hills per replicate at heading and 
at 20 days after heading. Whole plants with the normal number of panicles were 
collected from 8 hills per replicate at heading and at maturity. Samples were dried in 
an oven at 80 J C for 5 days, weighed, and powdered in a mill (W-100; Ikeda Rika, 
Tokyo, Japan). The nitrogen concentration on a dry-weight basis in each sample was 
determined with the Sumigraph NCH-22 analyzer (Sumika Chemical Analysis 
Service Center, Tokyo, Japan), and nitrogen content was calculated from the dry 
weight and nitrogen concentration. 

Measurement of xylem exudation rate and cytokinins in xylem exudate and flag 
leaves. Five hills with the normal number of panicles were collected per replicate at 14 
days after heading in 2011, and flag leaves were collected from the main tillers. To 
obtain exudates, shoots were cut 10 cm above the soil surface, and a polyethylene bag 
containing cotton wool was attached to each cut end with tape. Bags were detached 
3 h later, sealed and weighed. The weight of the exudate was calculated by subtracting 
the weight of the bag and the cotton wool. Exudates were collected by squeezing the 
cotton wool using syringe and immediately stored at — 80°C. Cytokinins were 
extracted using an automatic liquid handling system (MICROLAB STAR, Hamilton, 
NV, USA) for solid phase extraction and quantified by ultra-performance liquid 
chromatography (UPLC) coupled with a tandem quadrupole mass spectrometer 
(qMS/MS) equipped with an electrospray interface (ESI; UPLC-ESI-qMS/MS, UPLC- 
Xevo TQ-S, Waters, MA, USA) as described previously 38 . Cytokinin flux was 
calculated from the exudation rate and cytokinin concentration. 

Measurement of the pushing resistance of hills. Pushing resistance can be used as an 
index of lodging resistance of rice plants 27,39-42 . A prostrate tester (DIK7400; Daiki 
Rika Kogyo Co., Saitama, Japan) was set perpendicularly to the shoot at 15 cm above 
the soil surface, and pushing resistance was measured by pushing the plants to an 
angle of 45 from the vertical. Eight hills (with fertilization only) with the normal 
number of panicles per replicate were selected at 14 and 35 days after heading in 2011. 
The average number of panicles was 13.5 hill -1 for IR64 and 13.6 hill -1 for Drol-NIL 
at 14 days after heading, and 13.8 hill -1 for IR64 and 13.7 hill -1 for Drol-NIL at 35 
days after heading. Water in the paddy field was drained starting from 25 days after 
heading. 

Statistical analysis. Statistical analyses were performed by using a general linear 
model in SPSS version 17.0 (IBM, Tokyo, Japan). Analysis of variance (ANOVA) was 
used to test the effects of line and fertilizer treatment on root traits, yield, yield 
components, dry weight of the above-ground parts, and parameters of nitrogen and 
cytokinin. Line and fertilizer treatment were considered as fixed effects. Year and 
replication were considered as random effects for the former four parameters. 
ANOVA was also used to test the effects of line on pushing resistance. Replications 
were considered as random effects. Significant treatment effects (P <0.05 or 0.01) 
were explored by using Fisher's protected least significant difference test. 
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